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Combined static and molecular dynamics first-principles calculations are used to identify a direct 
structural link between the met ast able crystalline and amorphous phases of Ge2Sb2Tes. We find 
that the phase transition is driven by the displacement of Ge atoms along the rocksalt [111] direction 
from the stable-octahedron to high-energy-unstable tetrahedron sites close to the intrinsic vacancy 
regions, which give rise to the formation of local 4-fold coordinated motifs. Our analyses suggest 
that the high figures of merit of Ge2Sb2Tes are achieved from the optimal combination of intrinsic 
vacancies provided by Sb2Te3 and the instability of the tetrahedron sites provided by GeTe. 
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Ternary (GeTe) m (Sb2Te3) n materials, in particular 
the Ge2Sb2Te5 (GST) composition, have been considered 
as the most natural candidates for non-volatile memory 
applications through exploiting the fast and reversible 
resistance change between a metastable (m-GST) crys- 
talline phase (low resistivity) and an amorphous (a-GST) 
phase (high resistivity). [2,0, 0, 5[ However, the mecha- 
nism of the phase transition is still under intense debate. 
The existing models, [H|, 0, 0, Q have provided a prelimi- 
nary understanding of the transition mechanism, but fail 
to provide a clear and direct structural link between the 
m-GST and a-GST phases, which play a key role in the 
understanding of the reservible transition at an atomistic 
level. The m-GST g M M E3, M M phase crystal- 
lizes in a rocksalt-type (RS) structure, in which the Te 
atoms occupy the anion sites and Ge, Sb, and the nat- 
urally occurring intrinsic vacancies from Sb2Te3 (20% in 
GST) occupy the cation sites. It has been suggested that 
a-GST is characterized by the presence of 4- fold coordi- 
nated Ge atoms, [1, 0, [lfj E3, EH in which the sum 
of the occurrences GeTe4, Ge(SbTes), and Ge(GeTes) is 
about 66%. [19[ Ge— Ge and Ge— Sb bonds are found in 
those motifs, which is assumed to be due to disorder ef- 
fects, since they are not present in the crystalline phases. 
Furthermore, a-GST shows a volume ex pan sion of 6 — 7% 
compared with the m-GST phase, [10|, [20[ and it has a 
higher energy (28 — 40 meV/atom) with respect to m- 
GST.[2l[ Theoretically, first-principles molecular dynam- 
ics (MD) starting from a liquid phase with slow cooling 
rates have been used to generate a metastable phase (RS- 
type structure), however, no direct transition path was 
identified to link the proposed m-GST and a-GST phases. 
Thus, a new approach that connects the two phases at 
the atomistic level becomes highly desirable. 

In this work, using first-principles methods, we will ad- 
dress the following open questions: Is there a dominant 
structure link between both phases? What are the roles 
of GeTe and Sb2Te3 in GST? We will show in this Letter 
that combined static (zero temperature) and MD (high- 
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FIG. 1: (Color online) Structure models of the GST phases, 
(a) Metastable crystalline GST (m-GST). (b) m-GST (shift) 
structure, in which the Ge atoms occupy the 4-fold tetrahe- 
dron sites with lowest energy, (c) Amorphous GST obtained 
at zero temperature (am-GST) using modified m-GST struc- 
tures (m-GST with Ge shift), in which the tetrahedron Ge 
sites were initially occupied, (d) Amorphous GST obtained 
by high temperature molecular dynamics DFT calculations 
(a-GST). The Ge, Sb, and Te atoms are indicated in green, 
blue, and red, respectively. 



temperature) first-principles calculations can explain the 
phase transition mechanism between the m-GST and a- 
GST phases. Moreover, our study shows that generat- 
ing the amorphous phase from a known crystalline phase 
provides a better understanding of the structural rela- 
tionship between both phases. Thus, it provides a new 
avenue for further study of amorphous materials phase 
change transitions. 

Our static total energy and MD calculations are based 
on the all-electron projected augmented wave (PAW) 
method [H, Q and density functional theory (DFT) 
within the generalized gradient approximation (GGA- 
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PBE)0 as implemented in VASP. [H,^ To represent 
the metastable phase (RS-type structure), we employ a 
hexagonal (2x2x1) unit cell, in which the Te atoms are 
stacked along of the [0001] direction. The MD calcu- 
lations were performed employing cubic and hexagonal 
cells with 108 to 126 atoms. The total energies and equi- 
librium volumes for all structures in both crystalline and 
amorphous phases were obtained by full relaxation of the 
volume, shape, and atomic positions of the unit cell to 
minimize the quantum mechanical stresses and forces. 

To understand the phase transition, we first estab- 
lished the crystal structure of the m-GST phase. [2^ as 
shown in Fig. [TJ The obtained structure is consistent 
with experimental results and provides new insights into 
m-GST. [1, i, [Tl|, EH, EH In this layered-structure the or- 
dered intrinsic vacancies separate the building block units 
(GST), in which the Ge and Sb atoms are intermixed in 
planes. All Ge atoms are 6-fold coordinated in m-GST. 
However, it has been reported that up to one fifth of the 
Ge atoms are 4-fold coordinated with Te atoms in a-GST 
(GeTe4), while the remaining Ge atoms form 4- fold mo- 
tifs with combined Ge, Sb, and Te atoms, [l^ We notice 
that tetrahedral Ge atoms can be obtained by shifting 
the octahedral Ge atoms in m-GST along the hexagonal 
c direction, i.e., there are two tetrahedron sites for each 
Ge atom, Fig. [TJ In order to identify the lowest energy 
tetrahedron sites, we calculated the energetics for the 
occupation of each site by Ge atoms. The lowest energy 
sites are located in the intrinsic vacancy regions, while 
the highest energy sites are located in the center of the 
GST building blocks, i.e., there is a strong preference for 
the four-fold Ge atoms to be located in or near intrin- 
sic vacancy regions. Assuming that all the Ge atoms are 
shifted from their octahedron sites and occupy the low- 
est energy tetrahedron sites, we find that 50% of Ge will 
shift from the octahedra to tetrahedra along the [0001] 
direction, while the remaining 50% shift along of the op- 
posite direction. The m-GST (shift) structure in which 
all Ge atoms occupy the tetrahedral sites according to 
the distribution of intrinsic vacancies and energy barri- 
ers is shown in Fig. [TJd, which leads to the formation of 
Ge— Ge bonds. This configuration is highly unstable, and 
the system will relax without energy barrier to a lower 
energy phase (see am-GST structure in Fig. [Tfc). 

To provide a more direct structural link between the 
m-GST and a-GST phases, we first generated a-GST 
structures using first-principles MD simulations at high 
temperatures, T, using the same approach adopted in 
previous a-GST studies. [1, [13, El, Q Secondly, we gen- 
erated several amorphous structures from modified m- 
GST structures, in which a percentage of the Ge atoms 
(100%, 75%, 50%, 25%) are shifted to the tetrahedron 
sites from the lower energy octahedron sites (m-GST 
with shift Ge). The goal is to show that amorphous 
structures obtained in this way (am-GST in Fig. [T]) can 
preserve most of the structural features present in the a- 
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FIG. 2: Pair-correlation functions of various GST phases. 
Amorphous GST obtained by molecular dynamic calculations 
(a-GST, black lines). Amorphous GST obtained from occu- 
pation of tetrahedron sites in m-GST and complete relaxation 
(am-GST, red lines). Meta-stable GST phase (m-GST, blue 
lines, scaled by 0.50). 



GST generated by conventional MD calculations (a-GST 
in Fig. [Tj), and therefore provide a direct structural link 
and solid evidence to support the mechanisms that deter- 
mines the phase transition from m-GST to a-GST. The 
amorphous structures obtained by both approaches are 
shown in Fig. [H To quantify our analysis, we calculated 
the pair correlation (PC) functions, which are shown in 
Fig. [21 For the a-GST structures, the PC functions were 
averaged over five structures, while the PC functions of 
the am-GST structures were calculated for ten structures 
with different initial occupation of the Ge tetrahedron 
sites; the structure that provided the best agreement with 
the a-GST PC functions is shown in Fig. O 

Our PC function analysis shows that for all the am- 
GST structures, the one in which 50% of the Ge atoms 
are shifted from the octahedron to the tetrahedron sites 
along the hexagonal [0001] direction and the rest 50% 



TABLE I: Bond lengths (in A) of Ge 2 Sb 2 Te 5 (GST) in the 
amorphous and crystalline phases. 
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Te atomic position (Ge = 0,0,0) 

FIG. 3: Potential energy path for atomic displacements of Ge 
atoms along of the rocksalt (RS) [111] direction of GeTe. (a) 
Distorted RS structure, (b) Perfect RS structure, (c) Long 
Ge— Te bonds zincblende (ZB) structure, (d) Graphite like 
structure, (e) Perfect ZB structure, (f) RS-layer structure. 



moves along the [0001] direction reproduces almost all 
features present in the PC functions of a-GST, although 
some minor differences still exist. Furthermore, even mi- 
nor features are well-described by both structures, with 
the formation of Ge— Ge bonds and cavity regions, both 
of which have been identified as key characteristics of a- 
GST.[l,[l6|, Em We observe that am-GST structures 
in which less than half of the Ge atoms are moved to the 
tetrahedron sites do not yield PC functions similar to the 
a-GST structures, instead they show strong similarity to 
the PC function calculated for m-GST (see Fig.[2j). Fur- 
thermore, we observed that only the am-GST structures 
in which the Ge atoms initially occupy four-fold sites in 
or near the intrinsic vacancies lead to structure properties 
in good agreement with the calculated MD a-GST struc- 
tures. Thus, it suggests that the location of the intrinsic 
vacancies plays an important role in the phase transition, 
which can be explained by the lower energy barriers for 
Ge displacements close to intrinsic vacancy regions. For 
the lowest energy m-GST structures, in which the intrin- 
sic vacancies are ordered in a plane perpendicular to c. 
However, at high temperature or under non-equilibrium 
growth conditions the intrinsic vacancies may distribute 
more randomly among the cation sites, which is expected 
to play an important role in the pattern of shifted Ge 
atoms from their stable octahedra. 

Our predicted results are in good agreement with avail- 
able experimental data. For example, using the calcu- 
lated equilibrium volumes for both phases, we obtained 
a density of 5.89 g/cm 3 (m-GST) and 5.35 g/cm 3 (a- 
GST and am-GST), i.e. the amorphization gives rise 
to a volume expansion, which decreases the density by 
about 9.20%. The experimentally observed expansion is 
on the order of 6.4%. [10[ The volume expansion upon 



amorphization is a consequence of the Ge atoms moving 
to the lower coordination sites in the a-GST structures. 
Therefore, the smaller volume deformation observed in 
the experimental sample may indicate that the amor- 
phization process in not complete, [1, [l9[ i-e. not all the 
Ge atoms are moved away from their stable octahedron 
sites. 

Comparison of the total energies reveals that the a- 
GST structure is about 140 — 182 meV/atom higher in 
energy than the lowest energy m-GST structure, which 
corresponds to the energy limit between the fully amor- 
phized (100% shift of the Ge atoms) and the ordered m- 
GST structure. Differential scanning calorimetry mea- 
surements obtained 28 — 42 meV/atom.[21[ We found 
that the calculated energy differences decrease by about 
30 meV/atom if the intrinsic vacancies become disor- 
dered in m-GST. Furthermore, the energy difference 
could be much smaller (e.g. about 50 meV/atom) if only 
a fraction of the Ge atoms undergo site transitions. This 
again suggests that full scale amorphorization of GST or 
a complete ordering of Ge, Sb, and intrinsic vacancies in 
m-GST may not be typical in the GST phases. 

The averaged bond lengths calculated for both phases 
are summarized in Table U along with available experi- 
mental results. [H, [l6|, [HI The calculated bond lengths 
deviate by about 3 — 6% from the experimental results; 
however, most of the error is due to the use of GGA in our 
calculations, which systematically overestimates the lat- 
tice constants by about 3%. Furthermore, it is important 
to notice that the nearest-neighbor distances are spread 
over a large range of values, e.g. Ge— Te is from 2.67 to 
2.94 A and Sb-Te is from 2.86 to 3.23 A. We found a 
contraction in the averaged Ge— Te bond lengths in a- 
GST of up to 10% compared with m-GST, e.g. Ge-Te 
decreases from 2.87 - 3.24 A (m-GST) to 2.67 - 2.94 A 
(a-GST), while experimental measurements obtained a 
decrease of about 12%. Similar trends exist for Sb— Te. 

To understand the relaxation effects introduced by the 
shift of Ge atoms from octahedron to tetrahedron sites, 
we calculated the potential energy path along the RS 
[111] direction for GeTe as a function of Ge shift from 
octahedron (perfect RS) to the tetrahedron (zincblende) 
sites. The results are shown in Fig. 03 As expected, the 
distorted RS structure has the lowest energy (54 meV 
lower than the perfect RS structure), in which the dis- 
tortion is driven by Peierls-type level repulsion near the 
band edge. Unexpectedly, the zincblende (ZB) structure 
in which the Ge atoms occupy the tetrahedron sites with 
bond angles (Ge— Te— Ge) of 109.47°, is not a local min- 
imum as would be expected based on the general trends 
for binary semiconductors. In fact, we found that the ZB 
structure relaxes without energy barrier to the 'graphite- 
like' or to the 'long-Ge— Te' structures, which have lower 
energies than the high-symmetry ZB phase. Thus, Ge at 
ideal tetrahedral sites are intrinsically unstable in GeTe, 
which drives the Ge atoms at tetrahedral sites in GST 
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to move away and adopt a variety of lower symmetry 
coordination environments. 

The variety of coordination environments found in the 
GeTe energy surface is remarkable. From Ge site occu- 
pation of (0.25,0.25,0.25) to (0.40,0.40,0.40), three struc- 
tures have similar energies, i.e. 'long-Ge— Te', layered- 
ZB, and layered-RS. In 'long-Ge— Te', the Ge atoms form 
three short bonds (2.77 A) and one long bond (4.68 A) 
with the Te atoms. However, Ge is only three- fold co- 
ordinated in the layered structures with bond lengths of 
2.76 A, which is 2.90% (14.51%) smaller than the short 
(longer) Ge— Te bond lengths in the distorted RS struc- 
ture. As the layered GeTe structures are lower in energy 
than the graphite-like phase and only about 100 meV/f.u. 
higher than the distorted RS structure, it indicates a 
strong tendency of Ge atoms to form four-fold motifs 
with three short Ge— Te bonds (about 2.76 A) and bond 
angles of about 90°. Similar results, e.g. short bond 
lengths and average bond angles of about 90°, are ob- 
served by our calculations for a-GST, which is also con- 
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sistent with previous MD results for a-GSTj8|, [17, 
well as by experimental observations. IR fl9j Therefore, 
the inherent instability of Ge at the tetrahedral sites, 
low displacement energy, and unique coordination pref- 
erences of GeTe plays an important role in the formation 
of a-GST. 

In summary, using first-principles calculations, we ob- 
tained a direct structural link between the meta-stable 
and amorphous GST phases, as well as the role of the 
parent compounds. The Sb2Te3 provides intrinsic lattice 
vacancies, while GeTe contributes its RS-type structure 
in which Ge displacements along the RS [111] direction 
can be realized at low energy cost. The instability at 
the tetrahedral sites leads to the generation of disordered 
GST structures in which the Ge atoms are mostly four- 
fold coordinated with three short Ge— Te bond lengths. 
As the displacement has the lowest energy near intrin- 
sic vacancy sites, our analysis suggests that a high de- 
gree of amorphization can be achieved most easily when 
the system has a composition of (GeTe)2(Sb2Te3), i.e., is 
consistent with the observation that GST has the high- 
est figure of merit of all Ge— Sb— Te compounds. More- 
over, we show that generating amorphous materials di- 
rectly from its crystalline counterpart provides a better 
approach to understand these type of phase transitions 
present in phase change materials. 
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